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ABSTRACT 

We present results from an analysis of a Chandra observation of the cluster of galaxies 
Abell 3581. We discover the presence of a point-source in the central dominant galaxy which 
is coincident with the core of the radio source PKS 1404-267. The emission from the intr- 
acluster medium is analysed, both as seen in projection on the sky, and after correcting for 
projection effects, to determine the spatial distribution of gas temperature, density and metal- 
licity. We find that the cluster, despite hosting a moderately powerful radio source, shows a 
temperature decline to around 0.4 r max within the central 5 kpc. The cluster is notable for the 
low entropy within its core. We test and validate the XSPEC pro jet model for determining 
the intrinsic cluster gas properties. 
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1 INTRODUCTION 

Abell 3581 is a nearby (z ~ 0.0218), richness class 0, cluster of 
galaxies. The central dominant galaxy IC 4374 has an ex tensive 
optical emission-line filament sys tem iDanziger & FocardHll988l 
Johnstone, Fabia n & Nulsenlll987l) like that seen around the cen- 
tral galaxy in more massive clusters (e.g. the Perseus cluster, 
IConselice et all200lh . 

Abell 3581 was previously studied in X-rays, using data from 
the RO SAT and ASCA satellites by Ijohnstone. Fabian & TavloJ 
1 1998) who found that it had a cool intracluster medium (ICM) 
with kT ~ 2 keV, and a short radiative cooling time throughout its 
core. Those data admitted a classical cooling flow with a significant 
amount of intrinsic a bsorption. 

Recently, (e.g. If^teisonetalJ |200l|, Ijojmstone et alJ I2002L 



IPeterson et all2003llKaastra et all2004ISanders et all2004l it has 
been shown that clusters which have short central cooling times 
have much less gas at temperatures below about a third of the am- 
bient cluster temperature compared with that expected from sim- 
ple cooling. This inconsistency between theory and observations 
has led to the proposal of many different mechanisms to reheat 
the cooling gas. One popular subset of these heating models ap- 
peals to energy input from the active nucleus in the central dom- 
inant galaxy, whi ch is usually found to harbour a radio source 
lBurnsetalll997l) . 

The central dominant galaxy in Abell 3581 hosts the power- 
ful radio source PKS 1404-267. Since this is a cool cluster with 
a powerful radio source at its centre, it might be expected that if 
active galactic nuclei (AGN) are able to heat the ICM in clusters, 



this source might show the effects more clearly than richer clusters. 
However, although we will show that the ICM does not have sig- 
nificant amounts of gas below a temperature of kT ~ 0.8 keV, there 
seems to be no other clear evidence for heating of the ICM by the 
AGN. 

We note that Abell 3581 has a similar mean temperature and 
luminosity to the cluster Abell 1983 which was recently studied 
with XMM-Newton by IPratt & Arnaud J2003h . The properties of 
the inner core of Abell 3581 are however quite different from those 
of Abell 1983. 



2 OBSERVATIONS 

We present a Chandra observation (Obsid: 1650, Sequence num- 
ber: 800118) of the cool cluster Abell 3581 which was made on 
2001 June 7 with the S3 back-illuminated detector in the ACIS-S 
instrument. 

Data processing has been done using the CIAO package avail- 
able from the Chandra X-ray Centre while spectral fitting used the 
XSPEC package (Arnaud 1996). The events file was reprocessed to 
apply the most appropriate gain file to this observation, which was 
made with a focal plane temperature of -120C, and to remove peri- 
ods of bad background flaring. This resulted in a total of 7165s of 
good on-source exposure from the nominal 7260s observation. The 
PI column of the resulting events file was corrected for the temporal 
dependence of the gain by using the CORR.TGAIN package made 
available by Alexey Vikhlinin 1 . 

Analysis of the cluster emission used data processed for best 
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background discrimination from the 'very faint' data mode. Analy- 
sis of the point-source used standard 'faint mode' processing since 
the point-source counts are moderately affected by pileup. 

Throug hout this paper, we adopt a redshif t for Abell 3581 of 
z = 0.0218 Ijohnstone. Fabian & Tavloilll99 8) and a cosmology 
that assumes a flat universe with Hq = 70km s and £2„,=0.3. We 
use N Wright's web page 2 to calculate that the angular diameter 
distance to this source is 90.9 Mpc and its luminosity distance is 
94.9 Mpc. The linear scale is 441pc per arcsec. Wherever we quote 
uncertainties these are at the lcr level, unless explicitly stated oth- 
erwise. In general we have used the chi-square statistic to assess 
goodness of fit and tests for significance of additional parmeters, 
but have used the XSPEC implementation of Cash's C-statistic to 
determine par ameter uncertainties and confidence regions 3 . 

Recently B arnes & Nulsenl 120031) have carried out a detailed 
investigation of the fluctuations in the Galactic foreground absorp- 
tion equivalent Hydrogen column density in the direction of this 
clus ter. They reference two absolute measurements from the litera- 
ture: IStarketa?]ll992t). using a 2-degree beamsize find = 4.5 x 

cm~ 2 , while lHartmann & Burtori 1 19971) using a 35 arcmin 
beam found = 4. 1 x 10 20 cm~ 2 . In their investigation of fluctu- 
ations in the value of Nu towards this clu ster, using new da t a from 
the Australia Telescope Compact Array. [Barnes & Nulsenl (2003) 
find that at the 5o level fluctations larger than ±6.2 x 10 19 cm~ 2 
are ruled out. Throughout this paper we adopt therefore the value of 
Nu = 4. 1 x 10 20 cm~ 2 (from the more compact beam observation) 
as the value of the Galactic foreground absorption. 

In Fig. □ we show the Chandra X-ray image of the central 
~3.5x3.5 arcmin region of Abell 3581 in the 0.3-3.0 keV band. 
The X-ray events were binned to have 1 arcsec pixels before being 
adaptively smoothed. A clear detection of an X-ray emitting active 
nucleus is seen close to the centre of the cluster X-ray emission. 

The presence of the radio source PKS 1404-267 at the centre 
of this cluster has been known since the 1960s, but the Chandra 
image shows, for the first time, that there is an X-ray emitting point 
source associated with it. In Fig.|5|we show a close-up of the X-ray 
image (which has been unsharp masked) with a map of the radio 
source obtained from the National Radio Astronomy Observatory 4 
archive overlaid in contours. This 700s 20cm (1477 MHz) radio 
image was taken on 1987 September 21 with the VLA in the "A" 
configuration. The peak flux density is 396 mjy/beam with a syn- 
thesised beam of 3.34x 1.14 arcsec in position angle —27 degrees. 
Clear cavities in the X-ray emission are seen to the East and West 
of the nucleus, coincident with the lobes of the radio source. 
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Figure 1. Chandra X-ray image of Abell 3581 in the 0.3-3 keV band. The 
data have been binned to 1 arcsec pixels and adaptively smoothed. 
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3 THE X-RAY POINT SOURCE 

The X-ray point source is located approximately at 
RA(2000)=14:07:29.8, Dec(2000)=-27:01:04. We extracted 
a spectrum from a circular region of radius 1.8 arcsec centred on 
the point source. The major component of the background for 
this source is the surrounding cluster emission, so in this case, 
we extracted the background spectrum from an annular region 
surrounding the nucleus with radius between 1.97 and 3.94 arcsec. 
Spectra were rebinned to have a minimum of 20 counts per bin. 

2 http://www.astro.ucla.edu/-wri ght/Co smoCalc . html 

3 http : //xspec . gsfc.nasa.gov/docs/xanadu/xspec/manual/manual 

4 The National Radio Astronomy Observatory is operated by Associated 
Universities, Inc., under cooperative agreement with the National Science 
Foundation. 



Figure 2. Unsharp masked X-ray image of Abell 3581 (grey scale) with 
the 20cm VLA map of the PKS 1404-267 radio source overlaid. There are 
12 contours in the radio map spaced equally in logarithmic steps between 
0.7 mJy/beam and 390 mJy/beam. 

Redistribution matrices and ancilliary response functions were 
computed using calibration files from the Chandra calibration 
database version 2.26. This release corrects a previous one-channel 
offset error in the FITS embedded function file and incorporates 
a correction for the degradation in low-energy quantum efficiency 
due to the build-up of a contaminant on the filter wheel. 

We have modelled the spectrum of the point source as a power- 
law modified, initially, only by Galactic absorption. Preliminary fit 
hujesi tits were input to the WebPIMMS 5 count rate simulator which 

5 http : //heasarc . gsfc .nasa.gov/Tools/w3pimms .html 
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suggested that pileup was important at approximately the 8 per cent 
level. This is at a level where pileup is starting to become signifi- 
cant, so all subsequent fitting of the point source wa s done includ- 
ing the XSPEC implementation of the pileup model of lDavislfcOOll) . 

The best fit using the pileup model applied to a power-law with 
the Galactic absorption column density fixed at 4.1 x 10 20 cm~ 2 
gives a chi-square value of 11.2 for 20 degrees of freedom. Allow- 
ing the absorption to be freely fitted does not significantly decrease 
the chi-square value. 

In order to estimate the power-law index and its statistical un- 
certainties we use Cash's C-statistic <Cashlll979l implemented in 
the XSPEC cstat command) since there are many PI bins with 
the minimum number (20) of counts. The best fitting power-law 
model has a photon number index of 2.28^q'H, a 2-10 keV flux of 
2.16 x 10 -13 ergcm~ 2 s and a 2-10 keV (rest frame) luminosity 
corrected for Galactic absorption of 2.35 x 10 41 ergs -1 . The 0.1- 
10 keV luminosity is 1.08 x 10 42 ergs -1 . 

We can calculate the expected luminosity of the point source, 
assuming that it is a black hole accreting according to standard 
Bondi theory lBondi|[l952l) with a radiative efficiency of 10 per 
cent. The bolometric luminosity is given by: 

8.50 x 10 26 M bh « e 
^acc = ^5 erg s , 

where Mj,h is the mass of the black hole in solar masses, n e is the 
electron density in cm -3 , and T is the accreting gas temperature 

in keV. 

[ he work of lBettoni et alj l2003h gives the mass of the black 
hole in PKS 1404-267 as being in the range 8.5 - 9.5 x 10 s M® 
while the temperature and density of our innermost spatial bin cal- 
culated from the pro jet model in section 4.2 are kT = 0.81 keV 
and n e = 0.066 cm -3 respectively. These values yield a bolometric 
luminosity of 6.9 x 10 43 erg s _1 . If we assume that most of the ac- 
cretion luminosity comes out in the 0.1-10 keV band then the ratio 
of observed to expected luminosity is 0.016. The nucleus is under- 
luminous by a factor of ~ 64. This is much larger than any expected 
bolometric correction (< x 10). We note also that part of the 'point' 
source may actually be the innermost parts of the radio jets, at the 
nucleus itself. 

We note that appropriate values of temperature and density to 
put into the formula for the luminosity are those pertaining at the 
accretion radius: 

„ GM hh 

"acc ~ o i 

C 
L s 

where c s is the sound speed, or 

_ 3.705 x 10- n M b h 

"acc — — • 

For the innermost measured value of T in our data R acc = 0.04 kpc, 
which is much closer to the nucleus than our inner most bin can 
resolve (5 kpc). The work of |Pi Matteo et alj <2003l) on M87 has 
shown that going in from 5 kpc to the central 1 kpc region the intr- 
acluster medium temperature drops by about a factor of two while 
the density rises by about a factor of three. If similar trends occur 
in PKS 1404-267 they would increase the expected accretion lumi- 
nosity, and the discrepancy in the observed value, by more than a 
factor of 8. 

We note that underlumino us accretion is a lso seen in M87 
(in the Virgo c luster) bv iDi Matteo et"ail <2003l) . NGC 6166 (in 
A bell 2199) bvlDi Matteo et alj 1200 lh . and the Sombrero galaxy 
bv lPellegrini et alj 120031)" However, in these objects the luminos- 
ity deficit is rather larger (a factor of ~ 10 4 ) despite the black hole 




Figure 3. Map of the entire S3 chip showing the concentric regions from 
which spectra were extracted. North is to the top, East to the left. 



mass being very similar to the one consi dered here. Poss i ble so lu- 
tions to the discrepancy are discussed bv lDi Matteo et all 12003). 



4 CLUSTER EMISSION 

Since the surface brightness distribution of the cluster emission is 
quite circularly symmetric we have extracted spectra from a series 
of concentric circular annuli centred on the nuclear point source, 
but excluding counts from that point source, and one other at: 
RA(2000)=14:07:32.6, Dec(2000)=-27:00:41. Fig^shows the re- 
gions used overlaid on an image of the cluster. The radio source is 
contained almost entirely within the central 5 kpc radius region. 

It is not possible to use on-chip data from the same obser- 
vation as background since the cluster occupies the whole of the 
S3 chip. We therefore reprocessed the standard background field 
which most nearly matched the observation date and instrumental 
configuration to have the same gain file as the target observation. It 
was not possible to apply the time-dependent gain correction to the 
background events since the time of observation has been removed 
from the standard background fields. 

Background spectra were extracted from identical chip re- 
gions of the reprocessed background file by reprojecting the detec- 
tor coordinates on to the sky coordinates of the target observation. 

Redistribution matrices and ancilliary response files for these 
spatially extended regions were made using the CIAO 3.0.2 tools 
MKRMF and MKWARF. Background spectra were extracted from 
the same regions of the standard background fields as were used to 
extract the source spectra. 
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Figure 4. Absorption, temperature and metallicity profiles as seen in pro- 
jection on the sky. The dashed line shows the Galactic column density. 



4.1 Projected quantities 

Single temperature mekal iLiedahl et all 1 19951) plasma emission 
models, incorporating freely fitting Galactic absorption were fitted 
(separately) to the spectra extracted from each annulus. The absorp- 
tion, temperature, metallicity and normalizations are free parame- 
ters for each annulus. In Fig.[4]we show the absorption, temperature 
and metallicity profiles of the cluster gas derived from these fits. 

The absorption measured in the annuli beyond 40 kpc from 
the nucleus does seem to be slightly systematically low compared 
with the Galactic column density inferred from 21cm observations 
which is shown as a dashed line in Fig. [4] This may indicate the 
level of uncertainty in the model for the correction of the contami- 
nant that has caused the degradation of the low-energy response in 
the Chandra data. We also note that the absorption column density 
increases roughly linearly from 75 kpc in towards the centre of the 
cluster by ~ 7 x 10 20 cm -2 



Ijohnstone. Fabian & Tavloill 19981) found that the central 3 ar- 
cmin region from the ASCA data required a column density of 



9_j x 10 cm . Such a large value would be consistent with in- 
ner 40 kpc or 90 arcsec of the Chandra data. However, if the instru- 
mental absorption has been overcorrected in the current analysis 
then a greater degree of agreement betw een the ASCA and Chan- 
dra data is to be expected. The work of iBarnes & Nulsei] 120031) 
has shown that the Galactic foreground 21cm emission is constant 
to within 0.6 x 10 20 cm~ 2 (5cr upper limit on the fluctuation) in this 
region of the sky so it is unlikely that variations in Galactic absorp- 



tion can offer a viable explanation for the increased absorption. The 
excess absorption seems to be centred on the cluster which would 
argue for a cluster origin rather than a Galactic origin. The nucleus 
does not however appear to require any excess absorption (see sec- 
tion 3). 

Looking now at the temperature profile, it is clear that the cen- 
tral region of the cluster is cool, at kT ~ 1 .2 keV and that the tem- 
perature rises by nearly a factor of ~ 2 to reach a maximum at 
75 kpc from the nucleus. In the one region exterior to this the clus- 
ter temperature drops again to near 1.9 keV, although not with high 
significance. We note that due to the fact that the cluster is not cen- 
tred on the S3 chip the outermost region is only a partial annulus. 

The metallicity in the innermost region has a low value 
of ~ 0.2 times the sola r value (assuming the solar values of 
I Anders & Grevesseil 989). In the region between 5 and 90 kpc from 
the nucleus we find a much higher value of ~ 0.6 — 0.7 times the 
solar value. Further out the metallicity falls quickly to < 0.2 times 
the solar value again. 

In order to assess whether the fall in temperature beyond 
90 kpc is cluster wide, we have used an accumulative smoothing 
/ contour binning technique due to J. Sanders (to be described by 
Sanders et al. 2004, in preparation) to define spatial bins of given 
signal-to-noise ratio for spectral fitting. This allowed us to drop 
the assumption of circular symmetry. Briefly, a surface brightness 
map was produced from the raw counts using a smoothing algo- 
rithm which accumulated counts around each pixel until a signal- 
to-noise ratio of 15 was reached ("accumulative smoothing"). The 
binning method ("contour binning") then grew bins along direc- 
tions where the surface brightness is closest to the existing value, 
until the signal-to-noise ratio exceeded 30. Finally, bins were con- 
strained to have an edge length less than or equal to three times 
the circumference of a circle of the same area as the bin. Signal- 
to-noise ratios are calculated taking into account noise in the back- 
ground spectrum which is extracted from identical regions in chip 
coordinates from the standard background files. 

Spectra were extracted from these regions and binned to 20 
counts per spectral bin. New response matrices and ancilliary re- 
sponse functions were made corresponding to these regions which 
were then fit independently using Cash's C-statistic. The resulting 
temperature map is shown in Fig. [5] 

It is clear now that the south western part of the cluster is 
systematically hotter than the northern part. Since the outermost 
annulus comes only from the northern part of the cluster it is not 
surprising that the outermost annulus has a temperature below the 
one immediately interior to it since that partial annulus covers the 
hotter region. 

In Fig. [6] we show a close-up of the temperature map of the 
central region of the cluster with the radio intensity contours over- 
plotted. There is some suggestion that the radio source avoids the 
the coolest gas, which may be a direct indication that the radio 
source heats the gas. The precise locations of the boundaries of the 
individual regions selected for fitting are however somewhat sensi- 
tive to the smoothing and signal-to-noise parameters of the contour 
accretion technique, and are required to follow the surface bright- 
ness of the cluster emission. 



4.2 Correcting for projection effects 

We have fit the annular spectra accounting for projection of counts 
from the exterior regions into interior annuli using the pro jet 
model in XSPEC, and assuming spherical symmetry. The volumes 
of the shells sampled by each annulus are calculated from FITS 



The Galaxy Cluster Abell 3581 as seen by Chandra 




14ti07m403 14h07m30s 

EA (2000) 




Figure 5. X-ray temperature map of Abell 3581. 



-27: :30 ■ 



-27: 01 :□□ 



-27: 01 :30 — 



-27:02:00 ■ 




~l 1 1 1 1 

14h0 77n34s 14h0 77n32s 14h07rn30s 14h07rn28s 14h07rn26s 

RA (2000J 



Temperature (keV] 

Figure 6. X-ray temperature map of the central region of Abell 3581 with 
radio intensity contours overlaid. 



header items which describe the inner and outer boundaries of the 
annuli. We note that the cluster emission extends beyond the outer- 
most annulus of our analysis so that the emitting volume associated 
with the outermost annulus is too small. This causes too much pro- 
jection on to the next annulus in, so that the fit to that annulus un- 
derestimates its brightness. The resulting effect is minimal beyond 




Figure 7. Temperature and metallicity profiles for Abell 358 1 derived using 
the pro jet model. 



the outer two annuli because the amount of projection is relatively 
small due to the steep surface brightness profile; see section 4.3. 

The model which we have fitted consists of a single temper- 
ature mekal plasma emission model in which the temperature and 
metallicity at each radius are freely fit. A single value of the fore- 
ground absorption column density is also allowed to fit freely and 
obtains a value of 3.8 ± 0.7 x 10 2l) cm~ 2 , consistent with that ex- 
pected from the HI 21cm measurements. 

In Fig. Q we show the temperature and metallicity profiles of 
the cluster gas derived from these fits accounting for projection. 
The temperature climbs steadily from 0.8 keV within the central 
5 kpc to around 2.25 keV between 60-90 kpc. We note that the 
intrinsic temperature of the central bin is substantially cooler when 
taking into account the hotter gas seen in projection against it. In 
the outermost annulus the temperature drops back to 1.5 keV. Since 
this region is only a part of an annulus, to the north of the cluster 
(Fig|3), it is not clear whether this decline in temperature is due to 
substructure or in fact a coherent property of the cluster. 

The metallicity is consistent with 0.8 solar thoughout most 
of the cluster observed in our dataset. The very innermost region 
(within 5 kpc, but excluding the nuclear point source) and the re- 
gion beyond 90 kpc both have much lower values of metal abun- 
dance, at around 0. 1 times the solar value. 

The drop in metallicity in the central bin could be due either 
to an extended power-law component, or to the iron-bias effect 
iBuote & Fabianlll999h . The iron-bias effect results from fitting a 
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Table 1. Parameters, for the central spatial bin, of models fitted taking into 
account projection effects. Model 1 is pro jet * phabs (mekal) . Model 
2 is projet * phabs (mekal + po). Model 3 is projet * phabs 
(mekal + mekal). In models 2 and 3 the power-law and second mekal 
components are only present in the innermost spatial bin. Temperatures are 
given in keV, abundances are relative to the solar value and column densi- 
ties are in units of 10 20 cm~ 2 . "D of F" is the number of degrees of freedom 
in the fit. 
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plasma containing more than one distinct temperature component 
with a single temperature model. 

We have tested for the possibility of there being a spatially ex- 
tended power-law or a second mekal component by adding them 
(separately) to the baseline projet * phabs (mekal) model 
(Model 1), but only in the central bin. The results of these fits are 
compared with the baseline model in Table Q For the test with 
the power-law component (Model 2), both index and normaliza- 
tion were left as free parameters while for the test with a second 
mekal component (Model 2) we allowed a freely fitting tempera- 
ture and normalization but tied the metallicity to that of the first 
mekal component. For both tests we found that the metallicity was 
very poorly constrained and preferred to fit to unphysically large 
values; we have therefore imposed an upper limit of twice the solar 
value on this parameter. Indeed the best fit value of the metallic- 
ity hits this limit, but we note that allowing much higher values 
of metallicity does not reduce the value of the chi-square statistic 
significantly. Our lower limits on the metallicity presented inQcor- 
respond to an increase of 1.0 in the C-statistic from the value fitted 
with the metallicity at twice the solar value. 

For the model including the power-law component an F-test 
suggests that there is a 6 per cent chance of obtaining such a signif- 
icant drop in chi-square if the power-law component is not present. 
For the model including the second mekal component an F-test in- 
dicates that there is only a 0.7 per cent chance of obtaining such 
a significant drop in chi-square if the second mekal component is 
not present. 

We conclude that a two-component mekal model is a signifi- 
cantly better fit in the central bin than one with only a single mekal 
component and is preferred over one with an extended power-law 
component. However, in both of the two-component models the 
metallicity is sufficiently unconstrained that it becomes consistent 
with that seen in the second to seventh bins, counting outwards 
from the centre. 

We note that these data are also able to constrain the mass 
profile of Abell 3581. The results from this work will be presented 
as part of a larger sample of cluster mass profiles by Voigt & Fabian 
(2005), in preparation. 

4.3 Testing the projection correction 

We have run a series of tests to check that the pro j ct XSPEC model 
reproduces what is expected when given known synthetic data. We 



find that it generally gives good results. In Fig. [3] the results of one 
test which is based upon our dataset are shown. 

We start by using the temperature, density and metallicity pro- 
files appropriate for Abell 3581, accounting for projection effects, 
as determined by the pro j ct model. Since the outermost spatial re- 
gion shows a significant temperature drop we have only used data 
from the central seven spatial regions in this test. In order to con- 
struct a smooth synthetic cluster the density profile was fitted with 
the sum of two beta profiles, while the temperature profile was fit- 
ted with the functional form of lAllen et alj 1200 ll) . In both cases 
the outermost point (of the seven) was omitted from the fit (see 
discussion at the end of this section). 

We then constructed a spherically symmetric cluster assum- 
ing the ideal gas equation of state. There were 10,000 radial bins 
equally spaced out to a radius of 2500 kpc, at which point the clus- 
ter was truncated. 

Spectra appropriate for the physical conditions at each radius 
were then calculated using the mekal plasma code. The emission 
was integrated along the line of sight to give a projected spec- 
trum at each radius in the cluster. Next the spectra were accumu- 
lated radially, in the same annular rings as used in section 4.2, and 
scaled to the luminosity distance appropriate to Abell 3581 using 
our adopted cosmology. These spectra were then input to XSPEC, 
multiplied with a uniform Galactic absorption component, and used 
to generate synthetic pulse invariant (PI) spectra assuming the same 
responses and integration time as the Abell 3581 observation. The 
counts were perturbed following poisson statistics. 

The synthetic spectra were then fitted with the pro j ct model 
in XSPEC using the same techniques as in Section 4.2. In Fig.|8|we 
have overplotted the synthetic cluster profiles for temperature, den- 
sity and metallicity (dotted lines) with those derived from fitting 
the PI spectra with the mekal and projet models, allowing for 
Galactic absorption (points with error bars). It is clear that the fit- 
ted temperature and density points generally match very well with 
the 'real' synthetic profiles. There appear to be no systematic resid- 
uals in the temperature or metallicity profiles. In the density profile 
we note that the outermost point derived from the projet fits is 
very significantly above the cluster profile. This is because there are 
counts projected on to this region that cannot be taken into account 
correctly. The pro j ct model requires the projected outer radius for 
each annulus to be associated with that data set in order to calculate 
the volume and therefore the projection factors. For the outermost 
region the volume associated with the counts is wrong in the sense 
that it is too low. This leads to an overestimate of the surface bright- 
ness and density in the outermost annulus. Although this causes too 
much emission to be subtracted off the next inner annulus, it can be 
seen that the effect is very small and the errors due to this effect 
quickly damp out. We therefore have confidence in the projection- 
corrected results presented in section 4.2, apart from the density in 
the outermost annulus. 

4.4 The issue of entropy 

In Fig. [5] we show the temperature, electron density and entropy 
profiles deduced from the fits accounting for projection plotted 
against the scaled radius. Many authors choose to scale the radial 
coordinate by i?20() me radius within which the mean cluster den- 
sity is 20 times the mean den sity of the universe. R200 has been 
shown bv lCole & Lacevl 1 19961) to be proportional to the virial ra- 
dius, which is the appropriate scaling quantity in hierarchical self- 
simila r models of structu re formation. We have used equation 9 
of lSanderson et al] 120031) as an approximation to the virial radius 
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Figure 8. Simulated cluster (dotted lines) and parameters derived using the 
XSPEC proj ct model (crosses). 



for our scaling. Using r c i uster =2 keV (a value which is between 
that measured in t he two outer bins of our data, an d within the 
range measured by Jo hnstone. Fabian & Tavlorll 19981 Table 1) we 
find that the virial radius is l.HMpc. lSanderson et al]fc003h have 
shown that this equation can significantly overestimate ^200 partic- 
ularly in small systems with large temperature gradients. Since we 
have used the maximum temperature in the cluster (at least within 
the region covered by our data) rather than an emission-weighted 
average our i?20() ought to be reasonably accurate. 

In calculating the entrop y we have used the definition given 
bv lLlovd-Davies et alj l2000l) . in which the entropy is 

5= 7XT 2/3 . 

The outermost spatial bin at Q.\R y is subject to two problems. 
Firstly, as described in section 4.3 the proj ct model overestimates 
the density. This translates into too small a value of entropy in that 
bin. In order to avoid this problem we have estimated the elec- 
tron density by fitting the other points in the density profile with a 
double-beta model and extrapolated this to the position of the final 
spatial bin. This yields an electron density that is a factor of 1.66 
times lower than calculated by projct for that bin. In Fig. [9] the 
higher electron density is that determined directly from the proj ct 
model while the lower value is our corrected value. The entropy 
point (S = 1 20 ± 1 3 keV cm 2 ) that we plot at 0. IR V is derived from 
this corrected density. 

Secondly, the drop in temperature seen in this bin may not be 



representative of the cluster averaged profile since it is only a par- 
tial annulus and there may be azimuthally dependent temperature 
structure within the cluster. If the azimuthally averaged tempera- 
ture profile were to remain flat at the average of the values in the 
sixth and seventh annuli, counting outwards from the cluster cen- 
tre, it would have a value of 2. 1 keV (however, we note that forcing 
this value into the projct model would reduce the temperature 
of those inner bins). Values of kT = 1.8 — 2.2 keV were found by 
Ijohnstone. Fabian & Tavloil ll998l) from fitting various models to 
the ASCA data from the central 3 arcmin radius region (0.07_K V ) in 
this source. Using the value of 2. 1 keV for the temperature would 
increase the entropy at 0. IR V to 167 keV cm 2 . We consider this to 
be an absolute maximum value for this quantity. We also note that 
the virial radius increases as the square root of the cluster tempera- 
ture so it is not strongly affected by changes in the adopted cluster 
temperature. 

The entropy profile for Abell 3581 (as a function of radius 
scaled by the virial radius) lies well below the pro file for clus- 
ters w ith similar temperatures as published bv lLlovd-Davies et alJ 
(2000) (taking into account the Sk/i -1 ^ scaling for the different 
assumed values of Hq), This result stands, despite the presence of 
the relatively strong Parkes r adio source in Abell 3581. The en- 
tropy profile for Abell 1983 IPratt & Arnau3l2003l) . which has a 
similar temperature {kT = 2.1 keV), has entropy values which are 
~ 40 ke Vein 2 higher in the region covered by the Chandra data 
on Abell 3581 (the entropy scaling facto r (1 + z)~ 2 T expec ted 
from self-similar formation models, e.g. IPratt & Arnaudl2003l is 
the same within 1 per cent for these t wo objects). Many a uthors 
(see e.g. references i n iBinnev & Tabodll995l iFabian et all l200lL 
IChurazov et alj feoOl) have recently argued that the mechanical 
heating from an active galactic nucleus with jets can balance the 
cooling from high density gas in the cluster and thereby prevent 
the formation of a classical cooling flow. Such heating might be 
expected to increase the entropy of the clusters, particularly where 
the gas temperature is cool and the radio source is strong, but this 
is clearly not the case in Abell 3581. 

IPonman et alj 120031) discuss the entropy of intergalactic gas 
in elliptical galaxies, groups and clusters. Their fig. 4 shows the 
gas entropy at 0. l/?200 as a function of temperature for 66 virialized 
systems. Below about 4 keV there is a wide spread in this fiducial 
entropy. The value of entropy at 0.1R V in Abell 3581 calculated 
from the data at that position in the cluster (with the corrected value 
for the electron density) is 1 20 ± 1 3 ke V cm 2 , or 1 67 ± 1 8 ke V cm 2 , 
if we assume the high value of the temperature in the outermost an- 
nulus. By either of these measures, Abell 3581 has an entropy at 
0.1R V significantly below that for the mean of clust ers a t 2 keV 
~ 233 ± 30 keV cm 2 , as shown in fig. 5 of IPonman et alJ 120031) . 
Fig. 4 in that same work shows there is a large range in entropy 
for clusters at 2 keV. Our measured value for Abell 3581 would 
make it equal to their lowest entropy cluster (Abell 262) at 2 keV. 
Abell 3581 emphasizes the overall spread in properties of low tem- 
perature clusters. 

A power-law fit to the entropy as a function of radius (S = 
kr a ) gives a slope a = 0.87 ±0.03. Further including a constant 
offset in the entropy at zero radius does not produce a significant 
reduction in the chi-square statistic as determined using an F-test. 
We note that this is a much flatter slope than expected for pure 
shock heating of the intrac luster gas wh ich is expe cted to give a 
slope of a = 1.1 flbzzi & Normarl200lh . However. IPonman et alJ 
(2003) have shown, in their fig. 3, that such slopes are usually only 
seen beyond 0.2R y , which is outside the region covered by our data. 
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Figure 9. Temperature, electron density, and entropy profiles determined 
from the pro jet model. The radial coordinate has been scaled by the virial 
radius, l.HMpc. There are two points plotted in the electron density pro- 
file at 0.1i? v . The upper point is determined directly from the pro jet value 
while the lower point is determined from an extrapolation of a double-beta 
profile fit to the seven inner points. The entropy point plotted at 0.1i? v is 
calculated using the extrapolated fit value for the density. See text for de- 
tails. 



Table 2. Parameters of fits to cooling flow region. Model 1 is phabs 
(mekal). Model 2 is phabs (mekal + mekal). Model 3 is phabs 
(mekal + mkef low) with the lower temperature fixed at 0.08 keV. Model 
4 is phabs (mekal + mkef low) with the lower temperature allowed to fit 
freely. Temperatures are given in keV, abundances are relative to the solar 
value, the mass deposition rate M is given in solar masses per year and col- 
umn densities are in units of 10 20 cm~ 2 . "D of F" is the number of degrees 
of freedom in the fit. 
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Figure 10. Allowed values of mass deposition rate as a function of low 
temperature cutoff in the cooling gas. Contours are plotted at A C-statistic 
= 2.3, 4.61 and 9.21 which correspond to 68, 90 and 99 per cent confidence 
regions for two interesting parameters. 



The flattening of this slope in the inner regions of clusters may be 
due to the cooling out of the lower entropy gas. 

4.5 A classical cooling flow? 

We extracted a spectrum from the inner 3 arcmin (or 79 kpc) of 
the cluster (but excluding the point source) for comparison with 
the ASCA spectrum of Johnstone, Fabian & Taylor 1 1998). This re- 
gion is expected to contain most of the mass deposition if a classi- 
cal cooling flow is present Ijohnstone. Fabian & TavloJl998h . We 
have fitted a range of models to this region; the results are presented 
in Table|2| Initially we fitted a single temperature mekal model af- 
fected only by (freely fitting) Galactic absorption (Model 1). This 
is our baseline model for the region. 

Adding in a second mekal component (Model 2) in which the 
metallicity is linked to the first component reduces the chi-square 
value to 303.9 for 203 degrees of freedom and gives a significantly 
better fit (F-test significance^ — 1.9 x 10 ). 

Next we have substituted a cooling flow model for the sec- 
ond mekal component (Model 3) in which the upper temperature 
and metallicity of the cooling flow are tied to that of the mekal 



component, and the lowest temperature in the cooling flow model 
is below the Chandra band. Although this model is not a better fit 
than the two-temperature model in terms of the chi-square statis- 
tic, it only introduces one further fit parameter instead of two for 
Model 2 and an F-test shows the rediction in chi-square for this 
model to have a greater level of significance (1 — 6.4 x 10~ 6 ). The 
mass deposition rate fits to 8.1 ± 1.3 M, ; , yr~' . (Note that for com- 
parison wit h previous values of m ass deposition rate published by 
Ijohnstone. Fabian & Tavloil ll998t) this value should be doubled to 
account for the different assumed value of Hq.) 

Finally, allowing the lower temperature of the cooling flow 
to be a free parameter allows chi-square to drop further to 302.8. 
There is an F-test probability of 0. 15 of obtaining such a large drop 
in chi-square if the additional fit parameter were not required. The 
allowed mass deposition rate is now 12.7 + ^ M Q yr" 1 . In Fie.fTol 
we show contours of constant C-statistic in the mass depostion rate 
/ lower temperature plane. The three contours have A C-statistic 
= 2.3, 4.61 and 9.21 which correspond to 68, 90 and 99 per cent 
con fidence regions for two intere st ing par ameters. 

Johnstone, Fabian & Tavlor 1 19981) performed a surface 
brightness deprojection of the ROSAT HRI data. Within the central 
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3 arcmin (110 kpc in Hq = 50, qo = used in that paper), a mass 
deposition rate of 45 — 50 M Q yr~' was found. Scaling these values 
to our current cosmology would give a value of 22 — 25 Mq yr~' , 
around a factor of three higher than fitted in this section when the 
gas is allowed to cool to very low temperatures. Simila r results hav e 
been found for many other clusters by e.g . IPeterson et alj 1200 ll) . 
IPeterson et ail 12003). Kaas tra et alj J2004). a lower cutoff temper- 
ature to the cooling gas proposed. At a cutoff temperature of around 
~ 0.8 keV agreement with the deprojection mass deposition rates is 
achieved, albeit with a poor chi-square value. We note that 0.8 keV 
is the lowest temperature directly observed in the pro jet fits of 
section 4.2. 

The problem remains as to why the gas appears to cool by 
about a factor of ~ 2.5 yet not any further and does not accumulate 
at the lower temperature. 



thus dense, and has had few recent mergers. Cooling in its high 
density core may absorb all the kinetic energy produced by the 
present active nucleus so that little is available for the rest of the 
gas. d) The dense core may have been introduced into the cluster as 
a result of a merger, or e) The dense cool core may have caused the 
moderately powerful radio source by more effectively fuelling the 
central black hole. Radio source heating as inferred for the Perseus 
cluster (by the dissipation of sound waves; iFabian et all2003h may 
also apply, and scale, to Abell 3581. 

The core of Abell 3581 appears to be similar to that of the 
Perseus cluster, except that it is about an order of magnitude less 
luminous and has a temperature that is about one third of that in the 
Perseus cluster. Somehow the core has avoided the heating, or at 
least the effects of heating, common to most other low temperature 
clusters and groups. 



5 DISCUSSION 

So far we have considered heating in Abell 3581 only from the 
point of view of balancing the cooling from the ICM in the core 
of the cluster. However, there is a more global property of clus- 
ters that shows that heating has been important. The observed lu- 
minosi ty/temperature (L r /T x ) relation for clusters, where L x °c T x 
(e.g. lArnaud & Evrardlll999t) ra ther than T} as expected if only 
gravitational heating is involved <KaiseJl98d) . indicates that some 
form of extra heating has taken plac e. The le vel corresponds to 
about 2-3 keV per particle (e.g. lWu et all200d) . Low temperature 
clusters and groups show a large dispersion in the L x /T x relation 
I Helsdon & Ponman 2000), which may mean that although heating 
is widespread, it is not the same in all objects. One possible origin 
for the heating, which has been suggested by many authors, is the 
radio source which is often associated with the central galaxy. 

We have measured the unabsorbed X-ray luminosity in the 
central 60 kpc of both the Perseus cluster and Abell 3581 (ex- 
cluding the point sources) and find values of 1.5 x 10 44 ergs~' 
and 1.4 x 10 43 ergs -1 respectively in the 0.3-10 keV band, so the 
Perseus cluster is factor of ten brighter than PKS 1404-267 in this 
region. By comparison, the 1.4 GHz radio luminosity of the in- 
ner two radio lobes in the Perseus cluster is 2.0 x 10 24 WHz _1 
while the 1.4 GHz radio luminosity of the lobes in PKS 1404-267 
is 2.5 x 10 23 WHz -1 , showing that the distributed radio power in 
the Perseus cluster is also a factor of ten brighter than in PKS 1404- 
267. 

Abell 3581 is therefore particularly interesting as it is a 
low temperature cluster with a strong radio source where heat- 
ing se ems not to have been effective, even close to the active nu- 
cleus. Rovchowdhurv et al. (2004) have suggested that heating by 
the AGN should also affect the entropy profile of a cluster beyond 
the coolest parts. Even if the radio source has offset radiative cool- 
ing in the very innermost regions, it is clear that it has not had a 
strong influence on the gas entropy at OAR200, compared with sim- 
ilar low temperature clusters. 

Various explanations for this are possible, a) The radio source 
may be very young compared with the radiative cooling time of 
hundreds of Myrs, in which case radiative cooling would have dom- 
inated in the past. However, no rapidly cooling cluster has yet been 
found elsewhere, b) The major heating phase in clusters may have 
occurred 5 or more billion years ago, beyond z ~ 1 , when the cen- 
tral active galaxy or galaxies were much more important. This may 
have been a variable process which has left some clusters without 
significant heating, c) Perhaps Abell 3581 is particularly old and 
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